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ABSTRACT 

Gyrosynchrot ron radiation field j from midly relativistic 
electrons in a magnetoactive plasma a -e asymptotically calculated 
using the Green tensor and the Fourier transform; Lion. These 
fields consist of the two components which correspond to the 
ordinary and exti aordinary modes. Taking into account these 
fields, the emissivities and the absorption coefficients from 
an arbitrary distribution of electrons are calculated in order 
to discuss the intensity , spectrum and polarization ol gyro- 
synchrotron radiation. In general, the transfer of electro- 
magnetic energy takes place along a direction different from 
that of the wave normal since tne radiation fields have a non- 
vanishing component along the direction of the wave normal. A 
consideration is given on the problem of radiative transfer in 
relation to the Stokes parameters. 


♦NAS-NRC Associate with NASA 


INTRODUCTION 


The problem of gy rosynchrotron radiation radiation from 
midly relativistic electrons in a magnetoactive plasma occurs 
in astrophysics, radio astronomy and magnetospheric physics. 

The radiation from these electrons moving along helical 
trajectories has some unusual properties due to th? cispersion 
and anisotropy of the magnetoactive plasma as considered by 
many authors (e.g., Bunkin , 1957; Eidman , 1958; McKenzie , 1964; 
Liemohn , 1965; Mansfield , 1967; Pakh o mov et al., 1962, 1933; 
Sakurai and Ogawa , 1969; Fung , 1969a). They have studied the 
gyrosynchrotron radiation fiells from an electron moving in a 
helical trajectory in a magnetoactive plasma from various 
points of view. Eid man (1958) and Liemohn (1965), for example, 
have solved the problem in a cold and collisionless magnetoactive 
plasma by applying the well known Hamiltonian method of quantum 
electrodynamics. By introducing the Green tensor and a Fourier 
transformation, Bunkin (1957) and McKenzie (1964) have found 
different asymptotic solutions of the radiation fields in a 
cold and collisionless magnetoactive plasma. Similar method 
have been further developed by M ansfield (1967) using Fourier 
transformations. In the case of a collisionless thermal magneto- 
active plasma, Pakhomov , Aleksin and Stepanov (1962, 1963) have 
tried to calculate an asymptotic formula for the radiation fields. 
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Recently, Rama ty and Lingen f elte r (1967) and Ramaty (1968) 

have studied the influence of the ambient plasma on synchrotron 

and gyros vnchrotron radiation from energetic electrons in order 

to interpret the suppression of type IV radio bursts at low 

frequencies. In this case, they have assumed the ambient 

plasma to be isotropic and homogeneous since it appears that 

uu » a; and u;., and , where uu . oj and iw„ are the angular 

p H p H 'pH 

radiation, plasma and gyro-frequencies, respectively. At 
ii^esent, this suppression of synchrotron radiation from 
relativistic electrons has been applied to interpret the ob- 
servational spectra of solar and stellar radio emission (e.g., 
Ram aty , 1968; McCray , 1967; Fung, 1969a; Fun g and Yip , 1966; 
Scheuer , 1965) . 

The physical condition which are actually encountered 

near large sunspot groups, however, seems to be quite different 

from the simple case where u> » u; and uu . In fact, the type IV 

P « 

radio bursts are usually emitted by midly relativistic electrons 
which are moving within a plasma medium where the strength of 
the sunspot magnetic fields is very high, say 1 000 gauss (e.g., 
Kiepen heuer , 1953; Ellison , 1963; Bray and Doughead , 1965). 

We must, therefore, consider the effect of magnetic fields in 
dealing with the radiation characteristics of type IV radio 
bursts (Sakurai, 1964, 1965, 1970), and so we must solve the 
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gy rosy nc hro tron radiation fields from midly relativistic 
electrons moving in helical trajectories in a magnetoactive 
plasma which take account of the influence of angular gyro- 
frequency . 

From this point of view, we deal with the problem of 
gy rosynchrotron radiation from midly relativistic electrons 
in a magnetoactive plasma by means of the method similar to 
that of Bu nk i n (1957) and Mc Kenzie (1964). In our method, 
an asymptotic property of radiation fields is investigated 
by using a Green tensor and a Fourier transformation. In 
this paper, we assume that the medium is homogeneous, cold 
and collisionless and immersed in a static external magnetic 
field. The absorption and polarization of gyrosynchrotron radia 
tion in this medium are considered in relation to the problem 
of radiative transfer. 

FUNDAMENTAL EQUATIONS AND THE DERIVATION 
OF RADIATION FIELDS 

a ) Fundam en tal Equations 

Maxwell’s electromagnetic equations in a magnetoactive 
plasma are given, by using the complex dielectric tensor [K] 
(Stix, 1962), as follows: 

curl E - i - B (2-1) 

c 

curl B = - i — [K]*E + — j (2-2) 

c c 

div E? = 0 (2-3) 


4 



div [K ] • E - 4np 


(2 -4) 


D - [K]'E, B - H , (2-5) 

where E f H, D and B are the electric and magnetic fields 
and the electric and magnetic inductions, respectively, 
all of which are assumed to change with exj? (-lout) p 

and c are the electric current and charge density and the 
speed of light, respectively. 

A magnetoactive plasma is electromagnet ica 1 ly dispersive 
and anisotropic and may be characterized by a complex 
dielectric tensor Ck 3 . This tensor is a function of the 
angular propagation frequency and of the angie between the 
wave normal and the static external magnetic field in case 
of axial symmetry (S t ix , 1962). It is necessary for the 
following analysis to assume that this tensor is Hermitian 

[K] - Ck]*, (2-6) 

where the asterisk denotes the complex conjugate and trans- 
pose. Furthermore, in our treatment, it is assumed \.hat 

the magnetic susceptibility of the medium vanishes so that 

-» -* 

B - H. 

By the use if the electromagnetic potential, A and 0, 
the electromagnetic fields are expressed as 


E - i - A - grad 0, (2-7) 

c 

H - curl A. (2-8) 
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Thus, fcqs . (2-1) to (2-6) reduce to 


curl curl A + -- tKigrad 0 - — =• K]A - — j (2-9) 

c 2 c 

c 

- grad (i-KJ grad 0 ) - div — 1 • A - 4np (2-10) 

c 

Since the potentials are not uniquely determined by 
Eqs . (2-7) and (2-8), a Coulomb gauge condition may be im- 

posed on them (e.g., Jackson , 1962; Sakura i and O gawa , 1969): 

div [K] A - 0. (2-11) 

In this gauge, 0 is just the static potential of the source 
charge and so only denotes the corresponding static field 

~P r - 

D - LKJ'grad <P . Consequently, the radiation fields E and 
s t 

—p 

H are derivable from the vector potential A alone. Thus 
the vector potential A satisfies the following equation: 

V 2 J + 4 Lk]-A - - * 3 , (2-12) 

z c 

c 

which will be solved in the following discussion. 

If we make use of a Green Tensor t-G(r,r )J defined in 

e 

an infinite space with no boundary, that satisfies the 
equation 

2 

curl curl [G J - Lk]-lGj = 6 (r - r^) , 

c 


(2-13) 


the formal solution for the vector potential (r, t) is 
expiessed mathematically as 

A - — j [G (r , r ') 3 J (r,f) dr dt., (2-14) 

c - 

r ' , t ' 

“♦ -* 

where r and r' are the position vectors from the null point 
of the coordinate system to the point of observation and to 
the location of the radiating electrons, respectively, t and 
t' are the times with respect to the point of observation 
and the location of the radiating electrons, respectively. 

Taking into account that the Green tensor which satisfies 
Eq. (2-13) is a function of the difference in the coordinate 
r - r’ alone, the Fourier transformation of LG J is given by 

[G ] - — - — ~ J LG (k,uu)] exp{i k*(r - r')} dk',(2-15) 
(2h) J 

• • ^ 

where LG (k,uu)j is the Fourier component of [G j . 

The Dirac 6-function is also Fourier-transformed as 

6 (r - r') - — — « J exp{i k* (r - r')}dk (2-16) 

(2n) J 

By the use of Eqs . (2-15) and (2-16), Eq. (2-13) is reduced 
to 

2 

k (k tG]) - k 2 [G] + ~ CK] - CG] = - 

c 2 


1 


(2-13) 


When a static external magnetic JieJd lies along 
the Z-axis in the Cartesian coordinate system, the complex 
dielectric tensor is given by (Stix, 1962) 


where 


[sc] -r s 


i D 


0 


ID 

S 

0 


0 

0 

P 


S - | (R ♦ L), D - | (R - L) 


R - , X X 

1-Y* L 1-rY ’ 


P - 1-X 


(2-17) 


and 


2 2 

X - wp /u> , Y - u> /uu. 

H 


•y substituting Eqs. (2-17) into (2-13) \ we obtain 


where 


[G (k,u>) ] - - tQ t ]/Q q1 


(2-18) 


Q - \ P (k 2 - k 2 ) (k 2 - k 2 ) 

O 2 li ll , 1 11 11 , 2 

c 


(2-19a ) 


2 

u, 1 

2 

a , 2 


2P 


1 [2SP ^ - (S + P) k 2 

c 


+[(S-P) 2 k^-4PD 2 k 2 ^ + 4P 2 D 2 } 

c c 


(2-19b) 

i 
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and the elements of [Q ] are 

Q - -k 2 k 2 + (k 2 + k 2 - S \) (k 2 . p iS_ ) 

xx y z x y 2 x 2 

c c 

„ 2 2 

2 ju 2 ul; 

Q - k k k + (k k - l D =7) (k* - P =7 ) 

xy xyz xy 2 ' x 2 

c c 

2 2 oi^ 

Q - k k (k k -ID -7) + k k (k ♦ k - S ~) 
tcz yzxy 2' xzx y 2 

c 7 c 

2 2 2 2 

Q * k k k + (k k + 1 D (k - P 

yx xyz xy 2 x 0 

c c 

Q - (k 2 + k 2 - S \) (k 2 - P - k 2 k 2 
y V yz 2 ' 2 xz 


ID 




Q - k k (k k + 1 D — 7 ) + k k (k + k - S ~) 


yz 


zx 


zy 


zz 


x z x y 


c 

2 


y z 


c 

2 


m 2 2 

k k (k k + i D —7) + k k (k + k - S ~) 
yz xy 2 xz x z 2 

c c 

2 2 

k k (k k -ID ~) + k k (k 2 + k 2 - S 
xzxy 2 yz y z 2 

c c 

2 2 

, 2 2 <d % „ 2 2 uj v 

(k + k - S ~) (k^ + k - S ~) - 
y z 2 x z 2 

c c 

(k k -i D (k k + i D ~),k 2 - k 2 + k 2 + k 2 

\ y 2 x y 2 ’ xyz 

c c 

2 2 2 

and k x ™ k + k , k - k . 

x y 11 z 


Using Eqs. (2-15), (2-18) and (2-19), Eq. (2-14) is 
calculated as follows: 


A (r.t) 


. 3 
4rr c 


; 


c [Q^J J 7 (rlt") 

“ 2 “2 . 2". .,2 2 7 

P u> (k -k _ ) (k -k 0 ) 

11 u,l a 11 , 2 


exp [i (k.r-r') - u> (t-t')}] 
x dk du> dr' dt * 


( 2 - 20 ) 


This nives the vector potential due to the current density 


j (r ,t ). In order to calculate the gy rosynchrotron 
radiation from a mldly relativistic electron, we must define 
the current density by taking into account the helical 
orbit of this electron. 

b ) De r 1 vation of Gyrosynchrotron Radiation Field* 

The Cartesian coordinate system as shown in Fig. 

— ♦ 

1 is used here, where the static external magnetic field H 

o 

—p 

is along the Z-axis. The position vector r is assumed to 
be in the (y-z) plane and therefore r - (0, r sin or, 
r cos a) where cr is the angle between and the position 

-p -* 

vector. The wave vector k is decomposed as k - (k^ cos 0, 

—p — > 

k sin 0, k ), where 0 is the angle between H and k. 
x z o 

The orbit of a radiating electron in a helical trajectory 

—p 

is given as shown in Fig. 2, where r (t ) denotes the 

c 

position vector of the electron. Since the radiating electron 

—P -P t t 

produces the electric current, the current density j (r ,t ) 
in Eq. (2-20) is given by 


—* -p , . 

j (r » t ) 


e V 6 (r ' - r e (t 7 ) ) , 


(2-21) 


where e and V are the electronic charge and the velocity 
of the electron. The velocity V and the position vector 

—p . 

r (t ) of the electron with pitch angle f are expressed by 
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V 


(2-22a ) 


'JU uu 

- - .1 V, sin -r t ♦ a V cos — 7 t ♦ a V 

x x » y x T z n 


yV id yV uu 

r (t') - a — - cos — t' + a — - sin — t' + a' V t' ( (2-22b) 
e x u* H V y “h Y Z 


where V - V sin ♦, V - V cos ♦ , y - (1 

•L 11 


-»V* 


and fj - — . a , 
c x 


a and a are the unit vectors as shown In Fig. 2. 

y z 

By substituting Eqs. (2-21) into (2-20) and using Eq. 
(2-22), Eq. (2-20) can be integrated as follows ( Oga wa and 
Sakurai, 1969): 


A (r , t) - - 


e - c 2 [Q ]e*p{i(k A r A si n 0 + R^r 

— 5- J dkdui i-5 5- j 2 2 

2n c Ptt» (k -k ,) (k ,-k , ,) 

11 11 , 1 li 11 , 2 


- out)} 


£ [R„3 6 (iD-nuijj , - k u V u ) 


n»- 


(2-23) 


where 


[ H n 3 - -1 J n ^ (z)exp{ i(n-l) (0 - ^) } - J n+1 (z)expl i (n+l)(0 - ^) } 


v j. ,, p 

— (Jn_i(z )e x p {i(n-1) (0 - ) + J n + 1 (z) e *Pl 1 (n+l)(0 - 1 

V J (z) exp 1 n(0 - ^) 
u n -c 


k x V ^ 

and Z - Y . Here J (z) is the Bessel function of the first 

u> H n 
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kind of n-th order. The right hand side of Eq. (2-23) must 
be integrated with respect to k. In integrating this equation, 

we can decompose dk into k dk dk d 0. The integration over 0 
gives the radiation fields in the (y-z) plane. In order 
to obtain the radiation fields, all the terms of the order 
higher than r 2 are excluded in the integrand. The manipu- 
lative method for integrating over k is given in the 
literature (Ogawa and Sakurai , 1969). 

By integrating Eq. (2-23) with respect to k A , we obtain 
an asympototic expression of the radiation fields by applying 
the laborious method of steepest descents which will now be 
described. As has been defined in Fig. 1, k^ and k are 
given by 


k = k ( 0) sin 0 
J 


k - k. (0) cos 0, 
a J 


(2-24) 


where j 
cording 
waves in 
Ginzburg 


= 1,2 and 0 is the angle between H and k. Ac- 

o 

to the theory of the propagation of electromagnetic 

a magnetoactive plasma, k (9) is given as (e.g., 

3 


, 1964). 


k 2 (0) 

J 



2 

c 




4£ l e 3 


(2-25) 


where 


0 


2 2 
- P cos 0 + S sin 


e 2 - PS (1 + cos 2 6) + (S 2 - D 2 ) sin 2 0 


and 


c 3 " P (s2 " ° 2) 


Taking into account these relations between k. (j - 1, 2) 

%3 

and 0, we can transform the integration over k^ into one 
over 0. Thus, the integration over 0 leads to the following 
form: 


F (r) - J ♦ ( 0 ) exp{r$ (G)(d0 


(2-26) 


where c denotes the contour of the integration in the 

0 

0-plane and $ (0) is expressed as ik (0) cos (9-or). When 

the distance from the point of observation to that of the 

electron is much longer than the wave length of radiation 

fields, i.e., |k (0) r| » 1, the integration of Eq. (2-26) 

J 

is approximately carried out by the method of steepest 
descents. The result is 


F (r) 


L 


2rt 


r$" ( 0 ^ ) 


♦ (0 ) exp {ri (0 )} 


(2-27) 


where (6 ) 
s 


2 2 
b $/e 0 


and 0 is the angle at 
e s S 


the saddle point which can be determined by the equation 
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a. . t a k , (e) 

ae 


cos ( 6— or ) - ik (0) sin (0 - or) - 0, 

J 


from which we obtain the well-known equation, 


k . d 6 

J 


tan (9 - or) 


or 


i_ za 

. b 0 

J 


tan *, * - 0 - or 


which gives the relation between the wave normal direction 
0 and the ray direction o r (e.g., Stix , 1962; Bekef i , 1966), 
where (j » 1,2) defines the phase refractive index in the 
propagation of electromagnetic waves (^j - ck./w). 

Taking Eq. (2-27) into consideration, the asymptotic 
solution of Eq. (2-23) integrated with respect to k is given 


by 

A (r, t) 


i “V V (-1) J+1 j duu exp{i(— r cos ¥ - uut))l. (0) 
r * c j 

J n 


a .x 


u; 

^ H I 

X f- i U i 6 (uo-n — — — M’ u;0 cos 0)L _ 
• nj i Y J n 


V 

nj 


W 


nj 


f 2 - 28 ) 


where 
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I (0) 
J 


cos or(* sin 0) 

J 

COS X cos 0 


1 


[s in 



k ) cos X - 2 m- sin X) 
J J 


{ (S - 


P ) 2 M 


1 


sin 1 0 - 4PD 2 (a sin 2 0 - P)}]“^ 

• J 


U - 6 {^ 2 (P COS 2 0 + S Sin 2 6) - SP) J' (z) 

nj J- J n 


n “j. 

+ D { (P 

z 


* 2 Sin 2 6) 
J 


2 

3 m> sin 0 cos 0) J (z) , 
** J n 


V - 0 D (P - * 2 sin 2 0) J (z) 
nj a J n 


2 , 9 2 

+ (m - S) { (P - m sin 0) 

J z j 


- fj Mi sin 0 cos 0} J (z) 
li j n 


and 


W 


nj 


0 Dm sin 0 cos 0 J ' (z) 

x J n 


nH x 2 2 

- L m (m- . - S) sin 0 cos 0 

z J J 


+ P u {(j. 2 - S)(ii 2 cos 2 6 - S) - D 2 }] J n (z) 


Here 

M> . Y a) 3 sin 0 

* - -J — ,0 =0 sin *, 3 =3 cos ♦ , 

0) X 11 

H 

M - dp. /d0, M " " 6^M_/d6 2 and J* (z) = dj (z)/dz. 
J j J J n n 
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We must finally integrate Eq. (2-28) with respect 
to vu . Since the augument of the Dirac delta function in 
Eq. (2-28) is a function of the variable of integration, we 
therefore obtain the following relation, 


G(uu) 6 (f (uu) - a) ) duu 
b o 


G ( a) ) 
d f (uu) 

duu 


f (uu) 


(JU 

o 


(2-29) 


Since f(uu)*uu-uu^ 0 cos V cos 9 and uu 
L -®, »]) in our case, it follows that 

n V Y 

n 1 - £ cos V cos 0 


nuu / Y (n - 

n 


(2-30) 


which gives the well-known Doppler-shifted angular wave 
frequency . 

— > 

The vector potential A thus obtained is expressed as 


2 

A (r , t) - £ £ A (*.t) (2-31a) 

j = l n=-» J 


A . (r , t) 
nj 


'JU 

i ^ (-1)^ +1 A exp{i(— 


r cos ♦ 




V 

nj 


UU = (JJ 

n 



W 

nj 


(2 -3 lb) 


and 
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(2-31c) 


A 

j 


d (x M- ) 

1.(0) /| 1 - (3 COS 0 COS * — | 

J dx 


A (r,t) denotes the n-th 
nj 

vector potential. From Eq. 


tj. (x,0) is deduced to be 


-j _1 - 


2 2 

i y Y _ _s^ n 

I-X - 2 


harmonic of the j-th mode of the 
(2-25), the phase refractive index 


X U - X) 

4 4 

Y sin 


+ (1-X) 2 Y 2 cos 2 0)^ 


(2-32) 

where the upper sign (+), corresponding to the case j - 1, de- 
notes the ordinary mode of electromagnetic wave propagations 
and the lower sign (-) (j « 2), the extraordinary mode. 

As has been considered earlier in this paper, the radia- 
tion fields can be determined by using only the vector potential 
A. Thus the fields of the n-th harmonic of the j-th mode are 
easily calculated from the equations 


E-+1-A.B- curl A . (2-33) 

nj c nj nj nj 

c ) The Emissivi ties 

By calculating the Poynting vector, the flow of 
electromagnetic wave energy of the n-th harmonic of the j-th 
mode is given by 
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p -~(E xH*+E*xH ) 
nj 8n nj nj nj nj 


2 2 

e * 1*1 2 - 

cos X. A ( = n + ri a xn) 

2 J nj nj x id 


4ncr 


- ID 


(2-34) 


n 


e - e. 


s 


(j * 1,2) 


where the asterisk denotes the complex conjugate and 

2 2 

5 . - U . + z . 

nj nj nj 

2 2 2 2 
>; * (V - W ) sin cr cos cr + V . W (cos q - sin or.) 

nj nj nj J J nj nj J J 

Z * V cos a - W sin Of 
nj nj j nj j 

— * -+ 

n is the unit vector in the direction of r (Fig. 3). The 

first term of the Poynting vector is along the direction n 

which subtends the angle or from the external magnetic field 

J 

*+ . . 

as shown in Fig. 3, whereas the second term is in the (y - z) 

plane and perpendicular to the vector n. Since the second 
term does not give the energy flux passing through the point 
of observation along the vector r, we do not need to take 
into account this term in calculating the observational flux. 

The angular distribution of the radiation from a single 
electron at the n-th harmonic of the j-th mode at the point 
of the electron is thus calculated from Eq. (2-34) as follows: 
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Q nJ <e.0.*> 


2 2 

e ,j, x 



4rc 


o(>* x) 


cos x, A. { ,(1-3 cos 0 cos x 
J J nj 


ox 


J-> I 


X - X 


nj 


6 S (2-35) 


(j - 1,2) 


In deriving the above result, we have considered the correction 
discussed by Scheuer (1968) and Ginz bur g and Syrova tski i (1969). 
This result gives the emissivity from a radiating electron 
at the n-th harmonic of the j-th mode. Then, the frequency 
and angular distribution of the radiation from a single electron, 
into the ordinary and extraordinary modes, is given by 


r). (x,0,3,f) 


2 2 
e a; 

4nc 


d (m> x) 


cos X 


E A 5 . (1-3 cos 0 cos ♦ — ~ — ) 
“j J ro 3x 


(2-36) 


X - U) 


nj 


0-0 


s 

This result is very similar to that obtained by Ramaty (1969), 
but different from his case because the present result includes 
the effect due to the longitudinal component of radiation 
electric fields. This difference arises from the fact that 
the radiation energy does not propagate along the direction of 
the wave normal k because of the existence of the longitudinal 
component . 
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The emission-coefficient , j (w,0) are thus give*. by 

J 

the following general expression (Beke f i , 1966; Rama ty , 1969): 


J (uu , 0) - J T) (uu, 0,0,0 f ( p * ) d p’ ergs/sec str H cm' 
J J z 


(j - 1,2) 


(2-37) 


where f (p) d‘ p is the number of electrons per unit volume, 

3 -» 

with Hiomenta in d p around P. 

ABSORPTION 

The absorption coefficient, cr for an arbitrary 
velocity and pitch-angle distribution can be obtained from 
the following general expression (Bekefi, 1966): 


a (uu, 9,0 , t) 


3 2 

-|^ — J r\ (w, 0,0,O[f(p') - f(p)l d 3 p’ (3-1) 

ri -h u) J 

rj 


where h and r; are, respectively, Planck’s constant divided by 
2tt and the ”ray" refractive index given by Bekefi (1966). In 
the above equation, the sign is defined such that a photon is 
emitted when the momentum of an electron changes from p’ to p. 
Since the photon energy h<u is assumed to be much smaller than 


2 4 2 2 

the electron energy E (= vm^ c + p c ) , it follows that 


/“*, x ^ , \ E -h <ju df 

f (F ) - f ( P> “ 2 di + T~2 


-h a) 


pc 


p c sin ♦ cos ♦ 


x [E - 2 n, pc cos 0 cos ♦} ^ 
J d* 


(3-2) 


The coefficient ol the second term on the right hand side 
of the above equation gives the change in the pitch-angle 
of the radiating electron resulting from the emission ol a 
photon of energy into the direction 0. 

In case of an isotropic distribution (df/d¥ - 0), a 
necessary condition for negative absorption is df/dp > 0 
as can be seen from Eqs. (3-1) and (3-2). For an anisotropic 
distribution, however, the negative absorption can occur 
even if di/op <0. A necessary condition for this is given 
by 

— — * — {E - 2n . pc cos 0 cos ♦} -- > 0. (3-3) 

sin ¥ cos * 1 j * dl 


Since this equation contains the phase refractive index, n , 

« 1 

the condition for negative absorption is to some extent 
different between the ordinary ar.« extraordinary modes. For 
ultra-relativisitic electrons, the emission is strongly con- 
centrated into the direction of the instantaneous velocity and 


condition 

can be approx 

ima ted 



1 - 2n 

2 

cos 

♦ * f 



J 


df 



sin ♦ 

cos 

♦ 

and 

0 ~ ■ 

In the 

case 

of ♦ 

the 

pitch- 

-angle of 

radiating 


> 0 


21 


vtithin the range given by 0 < cos i < Vl/(2n ), in order to 
amplily synchrotron radiation. We must, however, remark that 
maser action from a moderately anisotropic distribution of 
radiating electrons will be confined mainly to the first 
few harmonics of gy rosynchrotron radiation (Hey vaer te , 1968). 

The momentum and pitch-angle distribution of radiatirg 
electrons is, therefore, very important in calculating the 
absorption coefficient and hence the source function in relation 
to the equation of radiative transfer of gy rosynchrotron radia- 
tion from midly relativistic electrons in a magnetoactive plasma. 

POLARIZATION AND RADIATIVE TRANSFER 
The shape of polarization ellipse, R n , for the n-th 
harmonic of the j-th mode is given by 



(k x a ) * E . (k x a ) * A 
x nj x nj 

— — ■ -*■ -» -> . 

a *A • k 
x nj 


a • E . * k 
x nj 


V . cos 0 - W . sin 0 
i — ^ ^ — — (J - 1,2) 


U 


nj 


- i (j - 1,2) 

nj 


(4-la) 

(4-lb) 


Taking into account that U and Z . are both real, the polariza- 

nj nj 

tion of both ordinary and extraordinary waves is elliptic with 
the axes of a and k x a / 1 k J (Fig. 4). 

X X 
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In a magnetoactive plasma, the electric fields E . which 

nj 

are given by Eq. (2-33) generally have a component along the 
direction k. In order to consider the radiative transfer in 
such a medium, we cannot, therelore, use the Stokes parameters 
defined for and applicable to the transverse polarized electro- 
magnetic waves. It is convenient to consider the electric 

-* — * 

induction D instead of using the electric field E in this case 

because the former does not have a component along the direction 

k (e.g., Zhe leznyakov , 1968; Fung , 1969b). After some algebra, 

— * 

the electric induction D is calculated as follows: 


D 

nj 


> , D , D , , ) 

nj x nj 0 nk k 


) (E 


_ F. 


r\ \ t j. __ i o\ 


(4-2) 


Here D , and E A are the components along the unit vector 
nj 0 nj 0 

a^ as shown in Fig. 4. The definition of the polarization, 

0 

Eq. (4-1), is thus rewritten as 


R 

nj 


E nj 0 
E 

nj x 


D nj 0 
D 

nj x 


Z 

= 1 U~ (J = lf2) (4_3) 

nj 


Since D does not have the component D . , , we can define the 
nj nj k 

Stokes p".draeters such as I , Q, U and V in an analogous way 


as the vacuum case, in which E . , =0 since V . sin 0 + 

nj k nj 

W . cos 0=0. 
nj 
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By adding both components of the ordinary and extra- 
ordinary modes, we can define the polarization tensor 
I ( o , 3 - x,0) of the radiation in terms of the components 

—p 

of D, which are given by 


D 

x 



X, 


D 


e 


E d . e 

nj 

j 


The tensor I _ is, therefore, defined as 
Cr A 


l a * - D a v <•. 


(4-4) 


where the asterisk indicates the complex conjugate of the 
components. The Stokes parameters are then expressed in 
terms of I * as follows: 


1 ! xx + I ee 


q - i - i 

xx 0b 


(4-5) 


u " X ex + I xe 


i (I rt - I _) 
0x x 0 


From Eq. (4-3), we can define the polarization coefficients 
in terms of the components of the electric induction of the 


radiation fields as 
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i 


(4-6, 



D 

xj 



1 °V (J 


1 , 2 ) 


Taking into account the ’’apparent" radiation intensities of 

the ordinary and extraordinary modes, 1^ and calculated 

from the electric induction and the phase difference 6 between 

both modes, D and D n can be written as 
’ x 0 


, D Xv /2 tTv^ fr i -io J n 

( D > " ( r ) [I 1 “l + 6 V “ 2 ’ 

0 


(4-7) 


where the unit vectors are given by 


a 


j 


(1 + “ej 





(j - 1,2) 


Consequently, we can express the Stokes parameters, Eq. (4-5), 

in terms of I and or (J - 1,2). In these expressions, we 

J 0J 

do not need to consider the component along the direction k 
of the electric fields. Therefore, we can develop the problem 
of radiative transfer in a similar way as has been considered 
by Zheleznyakov , (1968) and Fung (1969b). 

The degree of polarization, p, the ellipticity, e, and 
the angle X between the major axis and the x-direction can be 
written as (Born and Wolf, 1964; Bekefi, 1966): 
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p 


2 2 2 
V Q + U + V 

I 




tan 3 ; sin 23 


tan 2X. 


V 

2 2 

VQ + U + V 

U 

Q 


(4-8) 


Taking into account Eq. (4-5), the degree of polarization p 
is expressed as 


P 



(4-9) 


This means that p is always less than the unity when both modes 
can escape lrom the source region. 

As has been discussed earlier in this paper, the electro- 
magnetic energy emitted from electrons is usually transported 

in the direction making the angle c* with the static external 

J 

magnetic field which is not generally coincident with the angle 
as shown in Fig. 3. The Stokes parameters defined in this 

j 

paper are very useful for dealing with the problem of radiative 
transfer, especially the polarization treatment, but are not 
so much helpful in the study of the transport of the electro- 
magnetic energy in a magnetoactive plasma because of the reason 
just mentioned. 

The transport of the electromagnetic energy in a magneto- 
active plasma has recently been considered by B ekef i (1966) 
and Enome (1969) . According to them, the equation of the 
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transfer of electromagnetic energy is given by 


♦ * * 

+ . ( ii_) . 2l . a il_ 

ds ' 2* u) bt V 2 ; 2 ) 2 

n . g i n , n , n 

rj rj rj rj 


(j - 1,2) (4-10) 


* 

where the scalar quantity I ^ is known as the specific intensity 

of radiation and usually not equal to 1^ as defined earlier. 

Here dS is the element of the length along the ray direction 
— ♦ 

S (Fig. 3). ai is the group velocity. 

^ • ) 

In dealing with the transfer of radiation, let us assume 

-> *♦ 

the wave normal vectors, k and k to be initially in the 

I I 

same direction, i.e., k ||k . Even in this case, the ray 
direction of the ordinary mode in general, does not coincide 

with that of the extraordinary mode, , i.e., t a 2* 
result can be proved by calculating the equation 

tan (6 - <*.) - -i (j - 1,2) 

J . O 0 

J 

and 0=0-=* 0 9 . Accordingly, the two directions of the 

X m 

$ ♦ 

intensities I and 1 which are obtained by solving Eq. (4-10), 

X ^ 

are generally not coincident with each other. 

* 

The case in which I . and n are independent of time, 

J rj 

Eq. (4-10) is reduced to 
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OS ^ 7 ) 

n 

rj 


A 


- a 


A 
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rj 


j 2 
J n 

rj 


(j - 1,2) 


(4-11) 


or, by dividing both sides with or , 


OT 


° (-J-) 

v 2' 


J n 


rj 


n 


2 -v 


(j - 1,2) 


(4-12) 


rj 


where dt ( * - or, dS and 
j J 

* 

j I 

S - —A- (j - 1,2) (4-13) 

J n . j 

rj J 

which is known as the source function and plays an important 
role in the analysis of the radiation from the medium (e.g., 
Chandrasekhar , 1960). The above equation (4-12) shows that 
the effective path length is generally different between the 
two modes . 

When the phase difference between the ordinary and extra- 
ordinary modes is completely random, i.e., when the Faraday 
rotation is very large, the two modes of waves propagate 
independently. In this case, we can separately solve the 
equation of radiative transfer, Eqs. (4-11) or (4 -12) for 
each mode. For a homogeneous source region, we can integrate 
Eq. (4-11) as given by 
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(4-14) 


I 

j 


[1 - exp (- or L ) ], (j 

4 «J 


1 . 2 ) 


where L is the depth of the source. The intensity of the 
J 

escaping radiation can thus be given by 


"rj "j" [ ^ 1 (P,) dV 
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p sin v cos i 
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rj J 
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E df 
P dp 


p sin ♦ cos i|i 


{2 n. pc cos 
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G cos 1 


- E} — ) 


.3 ’ 
d p 


(j - 1.2) 


(4-15) 


In arriving at this result, we have used Eqs. (2-37), (3-1) 

and (3-2) . 

If the inequality 


E df. 
P dp 


p sin ♦ cos ♦ 


{2 


n 


pc i os 0 cos ♦ 


_ E) — 


< 0 


(4-16) 


is satisfied in some frequency range, the radiation can be 
amplified and then is observed to be enhanced in the same 
frequency range. 

The source function S is given, by using Eqs. (2-37), 

.3 

(3-1) and (3-2), as follows: 
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/ t) f (p' ) d p’ 


j I* . E df 


J P ^ 2 


p sin ♦ cos ♦ 


d f 3 

(2n.pc cos 0 cos * - E}— )d p' 


(J - 1,2) 


(4-17) 


The source function S is sometimes conveniently written 


in the form 


S - 
J 


1 it 


(X 


n 


2 or 


r~ 2 k V (J " x - 2) 


rj 


8n c 


(4-18) 


where k is Boltzmann's constant and T is a quantity with 

J 

dimensions of temperature. This is usually referred to as 
the radiation temperature of the medium and is a function of 
frequency, the direction of the r*»y in the medium and the 
elementary emission processes. The temperature is also 
usually different between the two modes, and is very important 
in estimating the radiation processes associated with solar 
and galactic radio emissions. This temperature occasionally 
becomes negative when maser action works in the medium where 


the radiation and propagation take place. 

I f we assume that u> » ix , ix and u> >:> u> , the refractive 

pH P H 

2 2 h 

index n is reduced approximately to (1 - tx /x ). If follows 
j P 

that n, - n„, n , = n 0 and x. - 0. - cr - 0. Taking into the 
1 2’ rl r2 j J j 

above assumption, the results obtained here are reduced to the 
same forms as currently used (e.g., Ramaty, 1969). It is clear 
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that the former are much more complicated than the latter. 
However, the equations of radiative transfer and their solution 
are similar to each other except for some essential differences 
on the solution of radiation fields. 

SUMMARY 

We have asymptotically calculated the gy rosynchrotron 
radiation fields from electrons in a magnetoactive plasma 
by usirg the Green tensor and the Fourier transformation. The 
medium has been assumed to be cold and collisionless and 
immersed in a static external magnetic field. The radiation 
fields obtained in this paper would be useful in studying 
the interpretation of solar radio type IV radio bursts and 
stellar and galactic radio spectra. 

Taking these radiation fields into consideration, we have 
developed the detailed calculation of the emissivity and 
absorption coefficients for a given direction of observation 
from an ensemble of electrons with arbitrary momentum and 
pitch an T le distribution. 

In order to study the problem of radiative transfer, we 
have used the electric induction D in place of the electric 
field E. The former is very convenient in formulating the 
Stokes parameters since it does not have a component along the 
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wave normal k. The electromagnetic energy is, however, 
usually transported in a direction different from the direction 
k. In order to calculate the transport of this energy, we 
must then consider the electric and magnetic fields, E and 
H. When the longitudinal component of the electric fields 
cannot be neglected, the theory of radiative transfer, therefore, 
becomes much complicated in comparison with the theory 
currently used. 

In this paper, we have not shown any results of numerical 
calculation, but this will be studied in a forthcoming paper 
and applied to problems on solar and galactic radio emissions. 
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